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I review the theoretical status of hadronization corrections to hadronic event ob- 
servables and discuss their impact on the determination of a s . 



Precise determination of the strong coupling has become one of the most 
important tasks of the QCD theory and phenomenology^. From the theoretical 
■ point of view the cleanest measurement comes from the total e + e~ annihilation 

cross section at high energies. The accuracy is dominated in this case by 
statistical errors; high statistics is needed because the effect proportional to 
a s is a few percent fraction only of the total cross section. Thus, going over 
to various event shape observables which are proportional to a s at the leading 
pH | order presents a clear experimental advantage. The price to pay is that the 

QCD description becomes more complicated and making a defensible estimate 
of the theoretical accuracy presents a nontrivial task. 

The uncertainty due to nonperturbative effects is a particularly delicate 



issue. From the experimentalist's point of view this is the uncertainty of 
'hadronization corrections' which are applied to uncover the structure of the 
event at the partem level from the observed structure at the hadron level. The 
hadronization process is modelled in a certain way, and differences between 
models (say, Lund string fragmentation or parton showers) are taken to esti- 
mate the error. In statistically important regions the hadronization corrections 
are of order 10% (at the Z peak) while the claimed error is of order 2-3%. This 
is less than the uncertainties of existing perturbative calculations (estimated by 
the scale dependence) so that for an experimentalist the hadronization could 
be considered under control. 

This procedure is very successful phenomenologically, but it is unsatis- 
factory from the theoretical point of view. In particular — I will return to 
this point later — the separation of perturbative versus nonperturbative alias 
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parton level versus hadronization effects is theoretically ill-defined, and the 
commonly accepted procedures might be suspected to be plagued by double 
counting of infrared effects. The theoretical understanding is guided by the 
Wilson operator product expansion (OPE) which is applicable, most notably, 
to the deep inelastic scattering (DIS). As a representative example, consider 
the Gross-Llewellyn Smith sum rule (GLS): 



Here ((O)) is the reduced matrix element of a certain quark-antiquark-gluon 
operator which quantifies the correlations between partons in the nucleon. 

From the OPE one learns: (i) Power counting of nonperturbative effects; 
in this case that perturbation theory is valid to 0(1 /Q 2 ) accuracy and that 
the description to order 0(1/Q 2 ) requires one dimensionful nonperturbative 
parameter; (ii) 'Universality' of nonperturbative effects in the sense that the 
same quark-antiquark-gluon operator appears in different physical processes; 
the coefficient is calculable in perturbation theory, hence one can sacrifice one 
measurement to get the prediction for other ones. 

Theoretical approaches to the nonperturbative corrections in hadron pro- 
duction would aim to get the similar structure. For a generic observable dom- 
inated by short distances one expects an expansion of the type 



The goal is to understand the power counting, that is which powers p arc 
present, calculate 'anomalous dimensions' 7 and relate the nonperturbative 
parameters Ar in different processes (universality) . Note that the question of 
actual magnitude of the nonperturbative parameters remains open. 

For the precision a s measurements the power counting is of primary im- 
portance. Indeed, assuming a 1/Q 2 correction, for the effective hard scale 
QcS > 10 GeV and the intrinsic size of 'matrix elements' 1 GeV, one obtains 
a ball-park estimate for the nonperturbative effects ~ 1% which can be ne- 
glected. On the other hand, if the nonperturbative effects are 0(1/Q), they 
can be of order 10% and one has to estimate them quantitatively. 

The structure of nonperturbative power-suppressed corrections to physical 
observables appears to be intimately connected with the same-sign factorial di- 
vergence r n ~ n! of the QCD perturbation theory in high orders. Divergence of 
perturbation theory implies that the sum of the series is only defined to power 
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accuracy exp[— const/a s (Q)] ~ 1/Q P and the ambiguity has to be remedied by 
adding nonperturbative corrections. In particular, the large-n behavior of the 
perturbative coefficients 

r n ~ const • n 7 • s~ n ■ n\, (3) 

necessarily requires a nonperturbative correction in (Q) with the powers p and 
7 related in a simple way to the coefficients s and 7, respectively. 

This has two consequences. First, one can investigate the structure of 
nonperturbative corrections to a large class of observables by studying the 
structure of higher orders_.in the perturbative expansions, which attracted a 
lot of recent activity (seecl for a recent review). Here is a short summary of 
results, related to a s determinations: 

• Most of the existing hadronic event observables are predicted to have 
nonperturbative corrections of order l/Qualx 

• The thrust and heavy jet mass distributions have l/Q corrections for 
the average values. However, outside the two-jet region nonperturbative 
effects have and extra suppression factor ~ a s (Q)U. 

• The one-particle inclusive cross section in the e + e~ annihilation (frag- 
mentation function) has only l/Q 2 corrections for fixed energy fraction; 
however, the integrated longitudinal (and transverse) cross sections have 
l/Q corrections LrO. 

• Power counting of the nonperturbative effects in jet fractions depends on 
the jet finding algorithm; l/Q corrections are intrinsic for.-JADE and are 
most likely absent for the Durham fc^-clustering method u. 

• There are no l/Q corrections to ihe Drell-Yan (heavy quark produc- 
tion) cross section to leading ordertd; their existence to the 0(a s (Q)/Q) 
accuracy is still disputed!!! 

The power counting of nonperturbative effects can be tested experimen- 
tally by the energy dependence of hadron event observables_subtracting the 
parton level prediction. This was done recently by DELPHIE3. 

Perturbative calculations in certain regions of phase space may require re- 
summations of large logarithms (threshold corrections). An important question 
is whether the power counting of nonperturbative effects is di st u rbed hy the 
resummations. It was studied for the Drell-Yan productionE30'E£rli3. Different 
resummation procedures which are equivalent in perturbation theory, can in- 
troduce different power-like corrections, and it was suggested that a criterium 



3 



for a 'good' resummation technique is that it does not bring in nonperturbative 
effects which are absent in finite .orders; this study emphasizes importance of 
large-angle soft gluon emission It turns out that quality of resummations 
in moment and momentum spaces is not related in any obvious way and their 
truncation may introduce spurious 1/Q P effects with a small power p ~ a s c3. 

The second consequence is that the separation of 'perturbative' and 'non- 
perturbative' effects is conceptually ambiguous. To make it meaningful one 
has to introduce an IR matching scale /^irEj, define 'nonperturbative' contri- 
butions as absorbing all effects from scales below ^ir, and subtract the small 
momentum contributions from perturbative series to avoid a double counting. 
Continuing the example with the GLS sum rule one obtains, schematically 

[ dxF^ p+pp (x,Q 2 ) = 3 
Jo 

(4) 

The premium for this refinement is that perturbation theory restricted to the 
contributions of scales above /Xtr is (almost) free from factorial divergences in 
high orders; the price to pay is that subtraction of small momenta is very awk- 
ward in practice and introduces an additional scale (and scheme) dependence. 

The same applies to the 'hadronization corrections'. The lessons to be 
learnt from OPE are: (i) separation of the parton cascade and hadronization 
is ambiguous; (ii) high orders of the perturbative series can imitate a power 
correction; (iii) if extracted from comparison with the data, the power correc- 
tion is expected to depend on the order of perturbation theory, factorization 
scheme and scale; (iv) numerical estimates suggest that 'true nonperturbative' 
contributions at small scales are of the same order as perturbative. 

To see how it works, consider energy dependence of the mean value of 
thrust (I — T) in the e + e~ annihilation. The experimental data are well de- 
scribed by the formula 

(I - T)(Q) = 0.335 a s (Q) + 1.02 o? s {Q) + (5) 

Here Q is the cm. energy, the two first terms on the r.h.s. correspond to the 
perturbation theory and the last term is the nonperturbative correction. The 
common procedure is to fix the 1 GeV/Q correction (assume it comes from a 
certain model) and fit a s (Q) to the data; this gives a s {Mz) — 0.120. Then, 
keeping the hadronization correction fixed, one varies the scale a s (Q) —> a s (/j,) 
to estimate the perturbative uncertainty coming from unknown higher orders. 

The caveat is that if nonperturbative and higher-order perturbative effects 
are inseparable, the hadronization correction can be scale-dependent itself, an 
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aspect which usually remains fogged. It should be fitted anew for each scale. 
Let us try \i = 0.13 Q which is of order of the gluon transverse momentum 
k± ~ (1 — T)Q 2 /A. Keeping a s (Mz) = 0.120 fixed, one obtains an equally 
good fit to the data by changing the 1/Q correction 

(1 - T)(Q) = 0.335 a,(0.13Q) + 0.19 a 2 (0.13Q) + (6) 

This illustrates that the reshuffling of higher-order perturbative corrections 
by changing the factorization scale can be compensated by the change of the 
hadronization correction. The only known way to separate the nonperturbative 
effects from the factorization scale dependence is to introduce an IR matching 
scale in the spirit of the OPE treatment as discussed above. This is attempted 
in the Dokshitser- Webber model where contributions of small scales are 
subtracted from perturbation theory. First applications of this model to the 
extraction of a s are encouragingl!3: The scale dependence is reduced by factor 
two compared to the traditional treatment, the IR matching scale dependence 
is small, and the value of the extracted nonperturbative parameter is stable. 

The question of possible 'universality' of hadronization corrections is com- 
plicated and needs further study. Real solution requires development of the 
OPE-likc techniques. The present approaches to this problem rely on two as- 
sumptions: that nonperturbative effects are proportional to perturbative am- 
biguities and that present 'naive' estimates of these ambiguities are representa- 
tive. This may well be not true, or only partially true, and these assumptions 
should be tested on simpler examples. Of particular interest are predictions for 
the x-deuendence of the 1/Q 2 effects in DIS □Ii3t3 and for the fragmentation 
functions LrB. The Dokshitser- Webber model has a buit-in universality also for 
the event shapes. 

In general, one should expect that the expansion is organized in powers of 
some physical scale (say the BLM scale) rather than the 'naive' cm. energy, 
and differences in the effective scales for various processes (together with the 
power counting) can give a rough idea of the difference in nonperturbative 
corrections. For example, the twist-4 effects in DIS at x — > 1 are proportional 
to A 2 /[(l — x)Q 2 ] while they are of order A 2 /[(l — x) 2 Q 2 ] for the Drell-Yan cross 
section. The difference reflects different hard scales for the gluon emission. 

Viewed this way, the problem of power corrections is inseparable from the 
familiar problem of the scale dependence and scale fixing. The specifics of 
hadron event observables is that the physical scale (estimated e.g. by any of 
existing scale fixing prescriptions) appears to be very low - of order 1/10 of the 
cm. energy. It was arguedca that the data do not show any preference for a 
low scale since the spread of the results for a s between different observables is 
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not reduced. From my point of view this conclusion is not warranted and the 
observed spread may indicate an intrinsic accuracy of the present treatment, 
with the separation of perturbative and hadronization effects as independent 
entities. {] 

In any case, some standartization is needed for the choice of scale in the 
data analysis. At present, each experimental collaboration uses different crite- 
ria to fix the scale range, which makes the comparison difficult. 

To summarize, what is done and what is necessary to do to improve the 
theoretical accuracy of the determinations of a s from hadronic event shapes? 

• Certain theoretical tools are developed to determine the power counting 
of nonperturbative effects in hadronic event observables. Some of the 
existing event shapes are predicted to have smaller corrections than the 
others, and it is advisable to concentrate on them rather than make 
'global fits'. 

• With this new knowledge, one should try to design new event shapes. 
An 'event shape of the year' should be measurable, calculable to O(al) 
accuracy and have small nonperturbative effects. Some work in this 
direction is reported in EL 

• In addition to the traditional procedures, one has to apply alternative 
methods for the data analysis with the hadronization corrections fitted 
to the data rather than taken from models. One should always bear 
in mind that hadronization corrections do not have objective meaning 
unless a clear scale separation is made. 

• The existing procedures should be checked for double counting of pertur- 
bative contributions at low scales. The complete next-to-leading order 
parton shower models would help a lot. 

• A convention for the scheme and scale setting is badly needed. 

My personal feeling is that present theoretical accuracy of a s (Mz) from hadronic 
event shapes is of order 7% and it can be improved up to factor two if the above 
questions are clarified. 

a Using low scales is more consistent with hadronization corrections taken from parton 
shower models which imply a s (k±) for each gluon emission. Note, however, that these 
corrections by construction parametrize contributions of small gluon virtualities, which is in 
spirit of the OPE. Combining them with the fixed-order perturbative calculations one has 
to subtract contributions of the low scales from the latter to avoid the double counting, at 
least in principle. In practice this is difficult to do because the parton showers do not include 
NLO radiative corrections consistently. 
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